Genomics in a changing arctic: critical questions await the molecular ecologist

Molecular ecology is poised to tackle a host of interesting questions in the coming years. The Arctic provides a unique and rapidly changing environment with a suite of emerging research needs that can be addressed through genetics and genomics. Here we highlight recent research on boreal and tundra ecosystems and put forth a series of questions related to plant and microbial responses to climate change that can benefit from technologies and analytical approaches contained within the molecular ecologist's toolbox. These questions include understanding (i) the mechanisms of plant acquisition and uptake of N in cold soils, (ii) how these processes are mediated by root traits, (iii) the role played by the plant microbiome in cycling C and nutrients within high‐latitude ecosystems and (iv) plant adaptation to extreme Arctic climates. We highlight how contributions can be made in these areas through studies that target model and nonmodel organisms and emphasize that the sequencing of the Populus and Salix genomes provides a valuable resource for scientific discoveries related to the plant microbiome and plant adaptation in the Arctic. Moreover, there exists an exciting role to play in model development, including incorporating genetic and evolutionary knowledge into ecosystem and Earth System Models. In this regard, the molecular ecologist provides a valuable perspective on plant genetics as a driver for community biodiversity, and how ecological and evolutionary forces govern community dynamics in a rapidly changing climate.     

Relocation,high‐latitude warming and host genetic identity shape the foliar fungal microbiome of poplars

Micro‐organisms associated with plants and animals affect host fitness, shape community structure and influence ecosystem properties. Climate change is expected to influence microbial communities, but their reactions are not well understood. Host‐associated micro‐organisms are influenced by the climate reactions of their hosts, which may undergo range shifts due to climatic niche tracking, or may be actively relocated to mitigate the effects of climate change. We used a common‐garden experiment and rDNA metabarcoding to examine the effect of host relocation and high‐latitude warming on the complex fungal endophytic microbiome associated with leaves of an ecologically dominant boreal forest tree (Populus balsamifera L.). We also considered the potential effects of poplar genetic identity in defining the reactions of the microbiome to the treatments. The relocation of hosts to the north increased the diversity of the microbiome and influenced its structure, with results indicating enemy release from plausible pathogens. High‐latitude warming decreased microbiome diversity in comparison with natural northern conditions. The warming also caused structural changes, which made the fungal communities distinct in comparison with both low‐latitude and high‐latitude natural communities, and increased the abundance of plausible pathogens. The reactions of the microbiome to relocation and warming were strongly dependent on host genetic identity. This suggests that climate change effects on host–microbiome systems may be mediated by the interaction of environmental factors and the population genetic processes of the hosts.     

Mapping Atmospheric Moisture Climatologies across the Conterminous United States

Spatial climate datasets of 1981–2010 long-term mean monthly average dew point and minimum and maximum vapor pressure deficit were developed for the conterminous United States at 30-arcsec (~800m) resolution. Interpolation of long-term averages (twelve monthly values per variable) was performed using PRISM (Parameter-elevation Relationships on Independent Slopes Model). Surface stations available for analysis numbered only 4,000 for dew point and 3,500 for vapor pressure deficit, compared to 16,000 for previously-developed grids of 1981–2010 long-term mean monthly minimum and maximum temperature. Therefore, a form of Climatologically-Aided Interpolation (CAI) was used, in which the 1981–2010 temperature grids were used as predictor grids. For each grid cell, PRISM calculated a local regression function between the interpolated climate variable and the predictor grid. Nearby stations entering the regression were assigned weights based on the physiographic similarity of the station to the grid cell that included the effects of distance, elevation, coastal proximity, vertical atmospheric layer, and topographic position. Interpolation uncertainties were estimated using cross-validation exercises. Given that CAI interpolation was used, a new method was developed to allow uncertainties in predictor grids to be accounted for in estimating the total interpolation error. Local land use/land cover properties had noticeable effects on the spatial patterns of atmospheric moisture content and deficit. An example of this was relatively high dew points and low vapor pressure deficits at stations located in or near irrigated fields. The new grids, in combination with existing temperature grids, enable the user to derive a full suite of atmospheric moisture variables, such as minimum and maximum relative humidity, vapor pressure, and dew point depression, with accompanying assumptions. All of these grids are available online at http://prism.oregonstate.edu, and include 800-m and 4-km resolution data, images, metadata, pedigree information, and station inventory files.     

AutoDockFR: Advances in Protein-Ligand Docking with Explicitly Specified Binding Site Flexibility

Automated docking of drug-like molecules into receptors is an essential tool in structure-based drug design. While modeling receptor flexibility is important for correctly predicting ligand binding, it still remains challenging. This work focuses on an approach in which receptor flexibility is modeled by explicitly specifying a set of receptor side-chains a-priori. The challenges of this approach include the: 1) exponential growth of the search space, demanding more efficient search methods; and 2) increased number of false positives, calling for scoring functions tailored for flexible receptor docking. We present AutoDockFR–AutoDock for Flexible Receptors (ADFR), a new docking engine based on the AutoDock4 scoring function, which addresses the aforementioned challenges with a new Genetic Algorithm (GA) and customized scoring function. We validate ADFR using the Astex Diverse Set, demonstrating an increase in efficiency and reliability of its GA over the one implemented in AutoDock4. We demonstrate greatly increased success rates when cross-docking ligands into apo receptors that require side-chain conformational changes for ligand binding. These cross-docking experiments are based on two datasets: 1) SEQ17 –a receptor diversity set containing 17 pairs of apo-holo structures; and 2) CDK2 –a ligand diversity set composed of one CDK2 apo structure and 52 known bound inhibitors. We show that, when cross-docking ligands into the apo conformation of the receptors with up to 14 flexible side-chains, ADFR reports more correctly cross-docked ligands than AutoDock Vina on both datasets with solutions found for 70.6% vs. 35.3% systems on SEQ17, and 76.9% vs. 61.5% on CDK2. ADFR also outperforms AutoDock Vina in number of top ranking solutions on both datasets. Furthermore, we show that correctly docked CDK2 complexes re-create on average 79.8% of all pairwise atomic interactions between the ligand and moving receptor atoms in the holo complexes. Finally, we show that down-weighting the receptor internal energy improves the ranking of correctly docked poses and that runtime for AutoDockFR scales linearly when side-chain flexibility is added.     

Dorsal Anterior Cingulate Thickness Is Related to Alexithymia in Childhood Trauma-Related PTSD

Alexithymia, or “no words for feelings”, is highly prevalent in samples with childhood maltreatment and posttraumatic stress disorder (PTSD). The dorsal anterior cingulate cortex (dACC) has been identified as a key region involved in alexithymia, early life trauma, and PTSD. Functional alterations in the dACC also have been associated with alexithymia in PTSD. This study examined whether dACC morphology is a neural correlate of alexithymia in child maltreatment-related PTSD. Sixteen adults with PTSD and a history of childhood sexual abuse, physical abuse, or exposure to domestic violence, and 24 healthy controls (HC) completed the Toronto Alexithymia Scale 20 (TAS–20) and underwent magnetic resonance imaging. Cortical thickness of the dACC was measured using FreeSurfer, and values were correlated with TAS–20 scores, controlling for sex and age, in both groups. Average TAS–20 score was significantly higher in the PTSD than the HC group. TAS–20 scores were significantly positively associated with dACC thickness only in the PTSD group. This association was strongest in the left hemisphere and for TAS–20 subscales that assess difficulty identifying and describing feelings. We found that increasing dACC gray matter thickness is a neural correlate of greater alexithymia in the context of PTSD with childhood maltreatment. While findings are correlational, they motivate further inquiry into the relationships between childhood adversity, emotional awareness and expression, and dACC morphologic development in trauma-related psychopathology.     

Longitudinal Stability in Reading Comprehension Is Largely Heritable from Grades 1 to 6

Reading comprehension is a foundational academic skill and significant attention has focused on reading development. This report is the first to examine the stability and change in genetic and environmental influences on reading comprehension across Grades 1 to 6. This developmental range is particularly important because it encompasses the timespan in which most children move from learning how to read to using reading for learning. Longitudinal simplex models were fitted separately for two independent twin samples (N = 706; N = 976). Results suggested that the shared environment contributed to variance in early but not later reading. Instead, stability in reading development was largely mediated by continuous genetic influences. Thus, although reading is clearly a learned skill and the environment remains important for reading development, individual differences in reading comprehension appear to be also influenced by a core of genetic stability that persists through the developmental course of reading.